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Abstract. In addition to scattering coefficients, the light ab- 
sorption coefficients of particulate and dissolved materials 
are the main factors determining the light propagation of 
the visible part of the spectrum and are, thus, important for 
developing ocean color algorithms. While these absorption 
properties have recently been documented by a few studies 
for the Arctic Ocean (e.g., Matsuoka et al., 2007, 2011; Ben 
Mustapha et al., 2012), the data sets used in the literature 
were sparse and individually insufficient to draw a general 
view of the basin-wide spatial and temporal variations in ab- 
sorption. To achieve such a task, we built a large absorp- 
tion database of the Arctic Ocean by pooling the majority of 
published data sets and merging new data sets. Our results 
show that the total nonwater absorption coefficients mea- 
sured in the eastern Arctic Ocean (EAO; Siberian side) are 
significantly higher than in the western Arctic Ocean (WAO; 
North American side). This higher absorption is explained 
by higher concentration of colored dissolved organic matter 
(CDOM) in watersheds on the Siberian side, which contains 
a large amount of dissolved organic carbon (DOC) compared 
to waters off North America. In contrast, the relationship be- 
tween the phytoplankton absorption ( )) and chlorophyll 


a (chi a) concentration in the EAO was not significantly 
different from that in the WAO. Because our semianalytical 
CDOM absorption algorithm is based on chi ^-specific a^iX) 
values (Matsuoka et al., 2013), this result indirectly suggests 
that CDOM absorption can be appropriately derived not only 
for the WAO but also for the EAO using ocean color data. 
Based on statistics, derived CDOM absorption values were 
reasonable compared to in situ measurements. By combining 
this algorithm with empirical DOC versus CDOM relation- 
ships, a semianalytical algorithm for estimating DOC con- 
centrations for river-influenced coastal waters of the Arctic 
Ocean is presented and applied to satellite ocean color data. 


1 Introduction 

The Arctic Ocean faces dramatic changes in physical en- 

vironments driven by ongoing global warming. Decreases 

in both sea ice extent and thickness (Comiso et al., 2008; 

Kwok, 2007; Stroeve et al., 2008) facilitate the propaga- 
tion of light into the water column, creating favorable con- 

ditions for phytoplankton’s photo synthetic processes when 
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nutrients are available in the upper layer of the water column 
(e.g., Ardyna et ah, 2013). Recent papers showed that an in- 
crease in primary production (PP) is significantly and pos- 
itively correlated with open water area (Arrigo and van Di- 
jken, 2011), suggesting that autochthonous dissolved organic 
carbon (DOC) resulting from the increased PP could be a 
nonnegligible source for open waters (Wheeler et al., 1996). 
More importantly, in terms of the DOC budget, DOC concen- 
trations in river-influenced coastal waters are much higher 
than for open waters (Benner et al., 2005; Stedmon et al., 
2011; Amon et al., 2012, Matsuoka et al., 2012). In addition 
to autotrophic sources, it is expected that recent increases in 
river discharges modify (Stedmon et al., 201 1) (1) the total 
amount of DOC discharged to the Arctic Ocean; (2) the sea- 
sonality of the riverine export by increasing the flux of DOC 
during the base-flow period (November-May), when mate- 
rial is older, more bio-refractory and could be transported far 
in the Arctic Ocean under the ice cover in winter; and (3) the 
relative contributions of the large Arctic rivers (e.g., Lena vs. 
Yenisei) with consequences on DOC distribution in the Arc- 
tic Ocean. However, a comprehensive method for quantify- 
ing and continuously monitoring DOC concentrations, which 
takes into account their temporal and geographical variabil- 
ity for river- influenced coastal waters of the Arctic Ocean, is 
presently not available. 

One of the best approaches to achieve this task is to ap- 
ply a semianalytical algorithm to satellite ocean color data. 
Such an algorithm was recently developed for the southern 
Beaufort Sea (Mackenzie Shelf area) waters in the Cana- 
dian Arctic (Matsuoka et al., 2013). The accuracy in the esti- 
mates of DOC concentrations using this approach relies upon 
an accurate retrieval of absorption coefficients of colored 
dissolved organic matter, CDOM, at 443 nm (<2cdom(443), 
m -1 ), and on the robustness of the relationship between 
DOC and <2cdom(443). The CDOM algorithm of Matsuoka 
et al. (2013) is based on a parameterization of absorption 
properties. To determine whether this algorithm can be ap- 
plied to regions of the Arctic Ocean other than the Macken- 
zie Shelf area, the geographical variability of those optical 
properties must be determined. 

The main objective of the present study is therefore to ex- 
amine absorption properties for various waters of the Arctic 
Ocean. Based on these results, a semianalytical algorithm for 
estimating DOC concentrations in river-influenced coastal 
waters of the Arctic Ocean is presented and applied to satel- 
lite ocean color data. 


2 Data sets and methods 
2.1 In situ measurements 

To provide a general overview of absorption properties, we 
built a large absorption database of the Arctic Ocean by pool- 
ing the majority of published data (Fig. 1 ; Table 1). Although 
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Figure 1. Locations of sampling stations for the SBl spr (green cir- 
cles), SBl sum (green diamonds), CASES (blue circles), MR aut 
(orange circles), NABOS (purple circles), MALINA (red circles), 
ICESCAPE1 (gold circles), and ICESCAPE2 (gold diamonds) in 
the Arctic Ocean. 


data from the Canadian Archipelago, the Greenland and Bar- 
ents seas were not included, our data set is quality checked 
and the largest for the Arctic Ocean. Additional data will be 
included to cover wider areas in the future. 

In the western Arctic Ocean (WAO), data were collected 
during the following five cruises from spring to autumn: the 
Western Arctic Shelf Basin Interaction, SBl, (1) spring and 
(2) summer cruises aboard the USCGS Healy (referred to 
as SBl spr: 5 May-15 June 2002, and SBl sum: 16 July- 
26 August 2002, respectively), (3) the Japanese Arctic cruise 
aboard R/V Mirai (referred to as MR: 1 September- 13 Oc- 
tober 2004), (4) the Canadian Arctic Shelf Exchange Study 
(CASES) cruise aboard CCGS Amundsen (referred to as 
CASES: 16 May-2 August 2004), and (5) the MALINA 
cruise aboard CCGS Amundsen (referred to as MALINA: 30 
July-27 August 2009). In the eastern Arctic Ocean (EAO), 
data were collected during the NABOS (Nansen and Amund- 
sen Basins Observational System) cruise aboard R/V Viktor 
Buynitsky (referred to as NABOS: 14 September-30 Septem- 
ber 2007). To minimize geographic bias due to a large num- 
ber of data in the Chukchi Sea and western part of the south- 
ern Beaufort Sea, absorption data sets from ICESC APE20 10 
and ICECAPE2011 cruises (referred to as ICESCAPE1: 15 
June-22 July 2010, and ICESC APE2: 25 June-29 July, 2011, 
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Table 1 . Summary of our absorption data sets. The total nonwater absorption (atw(A), m -1 ), and individual contribution by phytoplankton 
(<20 (A), m -1 ), NAP (tfNAP m_1 )> and CDOM GcdomW, m -1 ) are shown. To obtain a^iX) and ap(^) values, a different f 3 factor 
was used for each instrument because of the different geometry and phytoplankton species composition (see below; N indicates the number 
of samples). 


Year 

Season 

Cruise 

<7tw(443) 

(443) 

%a N AP(443) 

%<2cdom(443) 

N 

p factor 

2002 

Spring 

SBI spr 

0.10 ±0.04 

27 ± 0.15 

15 ± 10 

58 ±20 

68 

1 

2002 

Summer 

SBI sum 

0.11 ± 0.13 

31 ± 0.14 

16 ± 10 

53 ±16 

107 

1 

2003- 

2004 

Spring/Summer 

/Autumn 

CASES 

0.19 ±0.20 

16 ± 0.15 

15 ± 10 

68 ±18 

112 

2 a 

2004 

Autumn 

MR 

0.10 ± 0.06 

18 ± 12 

8 ± 5 

74 ±14 

179 

1 

2007 

Summer 

NABOS 

0.32 ± 0.15 

8 ±4 

8 ± 4 

85 ±7 

18 

3 a 

2009 

Summer 

MALINA 

0.19 ±0.45 

13 ± 9 

15 ± 15 

73 ±18 

119 

4 b 

Average 

All 

All 

0.14 ±0.23 

20 ±14 

0.13 ±11 

68 ±19 

603 

- 


1: Cleveland and Weidemann (1993) 

2: Tassan and Ferrari. (1995, 2002) 

3: Bricaud and Stramski. (1990); Allali et al. (1997) 

4: Rottgers and Gehnke. (2012) 

a Different factor was used for CASES (Tassan and Ferrari, 2002) and NABOS (Bricaud and Stramski, 1990; Allali et al., 1997) using the same PerkinElmer 19 
spectrophotometer equipped with a 60 mm integrating sphere. During the NABOS cruise (autumn), a mean value of chi b to chi a ratio of 0.15 was observed, 
which was similar to that obtained in the southern Beaufort Sea over the same time period (0. 19) when a mixture of small phytoplankton (e.g., flagellates) 
dominated (Matsuoka et al., 2009). This result suggests that small phytoplankton dominated in EAO waters during the NABOS cruise. Bricaud and 
Stramski. (1990) and Allali et al. (1997) obtained a ft factor for a smaller type of phytoplankton species, which was applied for the NABOS data set. For the 
CASES cruise, large phytoplankton such as diatoms dominated (Terrado et al., 2008). Tassan and Ferrari. (1995, 2002) determined a ft factor for larger 
phytoplankton, which was applied to the CASES data set. 

b A PerkinElmer 19 spectrophotometer equipped with a 150 mm integrating sphere was used. 


respectively) were not used for our main results. Instead, 
these independent data sets were used for evaluating the 
Bricaud and Stramski (1990) method to derive phytoplankton 
absorption (see Appendix A2). Radiometric data obtained 
from MALINA, ICESCAPE1, and ICESCAPE2 were used 
for evaluating a semianalytical CDOM algorithm developed 
by Matsuoka et al. (2013) (based on a method explained in 
Sect. 2.1.3 ). 

Discrete water samples were collected using Niskin bot- 
tles, except at the surface where a clean container was used 
during the MR cruise. During the six cruises, 30, 29, 51, 
94, 37, and 57 stations (SBI spr, SBI sum, MR, CASES, 
MALINA, and NABOS, respectively) were visited to collect 
water samples for pigment and absorption analyses within 
the euphotic zone (i.e., between the surface and the eu- 
phoric depth, z eu , defined as the depth of penetration of 
1 % of the surface photosynthetically available radiation 
(PAR)). Absorption coefficients of particles (a /? (A.), m -1 ) re- 
tained on GF/F (glass fiber) filters were determined from 
350 to 750 nm using a MPS-2400 spectrophotometer (Shi- 
mazu Corp.) for SBI spr, SBI sum, and MR, and a Lambda 
19 spectrophotometer (PerkinElmer) for CASES, NABOS, 
and MALINA. Because two spectrophotometers with differ- 
ent geometries were used for absorption measurements, the 
corresponding and appropriate pathlength amplification ef- 
fect caused by the highly scattering GF/F filter, known as the 
P factor (Mitchell and Kiefer, 1988), was taken into account 
for each instrument (see Table 1 for details). 

Phytoplankton pigments retained on the filters were ex- 
tracted using methanol (Kishino et al., 1985). Absorption co- 


efficients of nonalgal particles, NAP (<2 nap(^) ? m -1 ), were 
determined on the depigmented filter. Phytoplankton absorp- 
tion coefficients (^(A.), m -1 ) were finally obtained by sub- 
tracting &nap(A 0 from &p(X)- Fo r the NABOS cruise, only 
& P (X) data were available. To obtain a^iX) and ^nap(A-), a 
numerical decomposition method proposed by Bricaud and 
Stramski (1990) was applied to this data set. An evaluation 
of this method for Arctic samples is presented in Appendix 
A2. All absorption measurements were made following the 
NASA Ocean Optics Protocols (Mitchell et al., 2003) ex- 
cept the MALINA data, for which we followed Rottgers and 
Gehnke (2012). Spectral slopes of NAP, Snap (nm -1 ) were 
calculated by fitting a nonlinear model to the data from 380 
to 730 nm, excluding the 400^180 and 620-710 nm ranges to 
avoid any residual pigment absorption (Babin et al., 2003). 

For the CDOM absorption measurements, water samples 
were filtered using 0.2 pm pore-size filters immediately af- 
ter sampling (see details for sample preparation in Mitchell 
et al., 2003 and Matsuoka et al., 2012). Absorption coeffi- 
cients of CDOM (^cdom(A-), m -1 ) were determined from 
280 to 700 nm using a MPS2400 spectrophotometer (Shi- 
mazu corp.) for SBI spr, SBI sum, and MR. For the CASES 
cruise, the filtrates were frozen at —20 °C and shipped to the 
laboratory at Rimouski. (2 cdom(A-) spectra were then mea- 
sured using a Lambda 35 (PerkinElmer) from 250 to 800 nm. 
Both spectrophotometers were equipped with 10 cm quartz 
cells (Belanger et al., 2006; Matsuoka et al., 2011). For the 
MALINA cruise, an UltraPath liquid waveguide system was 
used (World Precision Instruments, Inc.), and ^cdom(A-) was 
obtained from 200 to 735 nm (Matsuoka et al, 2012). For the 
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NABOS cruise, ^cdom(^) was obtained at 412, 443, 488, 
510, 532, 555, 630, 676, and 715 nm using an ac-9 (WET 
Labs) with 0.2 pm inlet filter (Pegau et al., 2003). Spectral 
slopes of CDOM absorption (Scdom, nm -1 ) were calculated 
by fitting a nonlinear model to the data from 350 to 500 nm, 
except for the NABOS cruise, for which Scdom was calcu- 
lated using ^cdom(^) at 412, 440, 488, and 5 10 nm. The total 
nonwater absorption coefficient, a tw (X) was calculated as the 
sum of a<p( X), anapM, and acdomM- 

Chi a concentrations were determined using either fluoro- 
metric methods (chi a^ uo : Holm-Hansen et al., 1965; Suzuki 
and Ishimaru, 1990), high performance liquid chromatogra- 
phy (HPLC) (chi a HPLC : Ras et al., 2008), or both. A com- 
parison between chi < 2 fluo and chi a HPLC showed reason- 
able agreement (r 2 = 0.89, intercept = 0.11, slope = 0.98, 
N = 177; see Appendix Al). In this study, chi a HPLC was 
used if available. Otherwise, chi < 2 fluo was used. 

The detailed method for measuring concentrations of sus- 
pended particulate matter (SPM) is documented in Doxaran 
et al. (2012). Briefly, 0.2-6 L of seawater was filtered on a 
precombusted and preweighted GF/F filter. After filtration of 
seawater, filters were systematically rinsed with 0.2-0. 6 L of 
Milli-Q water to remove sea salt. The reason for rinsing using 
Milli-Q water instead of ammonium formate solution pro- 
posed by ICES (2004) is that the latter product may leave 
carbon deposits on the filters. This would have contaminated 
the CHN analyses, also made on the SPM filters to determine 
their particulate carbon content (see Doxaran et al., 2012). 
The concentration of SPM was then calculated as the differ- 
ence in weight before and after filtration divided by the vol- 
ume of the sample (gm -3 ). The SPM concentrations were 
measured for the CASES, NABOS, and MALINA cruises. 

2.1.1 Calculation of remote sensing reflectance 

In- water upwelled radiance (L u , pWcm -2 nm -1 sr -1 ) and 
downward irradiance (Ed, pW cm -2 nm -1 ) were obtained at 
19 wavelengths spanning from 320 to 780 nm following the 
NASA Ocean Optics Protocols (Mueller and Austin, 1995) 
and Hooker et al. (2013). Briefly, a compact-optical profiling 
system (C-OPS, Biospherical Instruments Inc.) (Morrow et 
al. 2010) was deployed at 36 and 19 stations for MALINA 
and ICESCAPE1 cruises, respectively. For the ICESCAPE2 
cruise, a profiling reflectance radiometer series 800 (PRR- 
800, Biospherical Instruments Inc.) was deployed at 24 sta- 
tions. In the present study, the independent data sets from 
these three cruises were used to evaluate a semianalytical al- 
gorithm for deriving CDOM absorption developed by Mat- 
suoka et al. (2013). The above-water global solar irradiance 
( E s ) measurements were used to correct the Ed and L u data 
for change in the incident light field during water column 
profiling. The in- water measurements were made far away 
from the main ship to minimize platform perturbations, and 
the above-water solar irradiance measurements were made 
to avoid superstructure shadows and reflections. The tilt an- 


gles were always less than 5°, as recommended by the NASA 
Ocean Optics Protocols (Mueller and Austin, 1995). Subsur- 
face L u at null depth (i.e., L u (0 - , X )) was obtained from 
the slope and intercept given by the least-squares linear re- 
gression of the log-transformed upwelled radiance versus z. 
The principal data product used here is the remote sensing 
reflectance, R r S (X) = 0.54 L u (0 - , X)/E s (X), where X indi- 
cates wavelength. In this study, R rs (X) at six wavelengths (i.e, 
412, 443, 490, 532, 555, and 670 nm corresponding approxi- 
mately to MODerate-Resolution Imaging Spectroradiometer, 
MODIS, ocean color bands) were used to derive <2 cdomW 
at 443 nm values based on a method explained in Sect. 2.1.3 
(Matsuoka et al., 2013). 

2.1.2 Ocean color data 

Level 3 monthly averaged climatology data over the 2002- 
2012 period for R rs ( X) at 412, 443, 488, 531, 555, and 
670 nm, obtained from the MODIS Aqua ocean color 
sensor (http://oceandata.sci.gsfc.nasa.gov/MODISA/Binned/ 
Monthly_Climatology/), were used in this study. A semi- 
analytical algorithm for deriving <2cdom(443) based on the 
method detailed in Sect. 2.1.3 (Matsuoka et al., 2013) was 
then applied to the R vs (X) data. DOC concentration in river- 
influenced coastal waters were estimated using the relation- 
ships between DOC and <2cdom(443) as shown in Sect. 2.1.3. 

Level 2 daily R rs ( X) at 412, 443, 488, 531, 555, and 
670 nm, obtained from the MODIS Aqua ocean color sen- 
sor (http ://oceancolor. gsfc .nasa. gov/egi/browse .pl?sen=am) , 
were used to generate daily DOC concentrations images and 
examine its spatial and temporal variability (see Fig. 13). 

2.1.3 Semianalytical algorithms for deriving CDOM 
absorption and DOC concentration 

A detailed description of our semianalytical CDOM absorp- 
tion algorithm is given in Matsuoka et al. (2013). Briefly, this 
algorithm produces <2cdom(443) using R rs (X ) at six wave- 
lengths in the visible spectral domain corresponding to ocean 
color sensors (e.g., MODIS and MEdium Resolution Imag- 
ing Spectrometer, MERIS) by minimizing the difference be- 
tween measured R rs (X ) and R rs ( X) calculated using absorp- 
tion and backscattering coefficients, as in the Garver-Siegel- 
Maritorena (GSM) algorithm (Garver and Siegel, 1997; Mar- 
itorena et a., 2002). The following modifications were in- 
cluded: (1) chi a specific absorption coefficients of phy- 
toplankton for Arctic waters (Matsuoka et al., 2011) were 
used; (2) spectral slopes of CDM absorption, Scdm (nm -1 ) 
and of backscattering of particles (b\> v (X), m -1 ), r) (dimen- 
sionless) were set either to be constants or as functions 
of R rs (X) (those parameters were determined based on a 
sensitivity analysis following Matsuoka et al. (2013); (3) 
Z?bp(555) was calculated using &bp(443) and (4) <2nap(443) 
was calculated using the empirical relationship between 
^nap(443) and bt p ( 555) obtained for Arctic waters (i.e., 
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<2nap(443) = Z?bp(555)/0.2393; Matsuoka et al., 2007); and 
(5) <2cdom(443) was obtained by subtracting &nap(443) from 
^cdom(443). To estimate DOC concentrations in the near- 
surface and river- influenced coastal waters using ocean color 
data, a relationship between DOC and <2 cdomW needs to 
be established. Walker et al. (2013) recently showed a con- 
sistent relationship between DOC and CDOM absorption 
coefficient at 350 nm (<2cdom(350), m -1 ) for waters at the 
mouth of the five major Arctic rivers (i.e., Lena, Yenisei, 
Ob, Mackenzie and Kolyma rivers; DOC (mgL -1 ) = 2.907 
+ 0.4059 • <2cdom(350)). The DOC vs. ^cdom(350) relation- 
ship was converted into a DOC vs. <2cdom(443) relation- 
ship by assuming a Scdom of 0.0175 nm -1 calculated us- 
ing data from 350 to 600 nm in the EAO, as reported by Aas 
et al. (2002) (i.e., DOC (pM) = 245 + 171 • ^ C dom(443)). 
Stedmon et al.(201 1) reported a lower Scdom (0.0167 nm -1 ) 
using the data from 300 to 650 nm. While the use of a differ- 
ent data range for calculating a spectral slope makes it diffi- 
cult to compare results (Twardowski et al., 2004), the spectral 
slope reported by Stedmon et al. (2011) would not be signifi- 
cantly different from the one by Aas et al. (2002) if the same 
data range is used: In any case, the choice of a spectral slope 
resulted in a change of only 7 % for the slope of the DOC 
versus «cdom(443) relationship, which did not influence our 
results. 

In the relationship converted in this study (DOC = 245 
+ 171 • £cdom(443)), the intercept is too high for the mouth 
of the Mackenzie River (245 pM) in the WAO, when com- 
pared to published values (55-97 pM; Osbum et al., 2009; 
Matsuoka et al., 2012). The high intercept may be influenced 
by a lower DOC to ^cdom(443) ratio at high &cdom(443) 
values (>l.lm -1 ) in the EAO (see also Fig. 10). Thus, 
for the WAO, we used the relationship recently obtained 
by Matsuoka et al. (2012) (i.e., DOC (pM) = 55 + 
357 • £cdom(443)). These two regressions were used to es- 
timate DOC concentrations from space using <2cdom(443) 
derived from ocean color data. Note that DOC concentra- 
tions were estimated only for river-influenced coastal waters 
where tight relationships between DOC and CDOM are ob- 
served. 





a lw (443)(m^ _1 ) a lw (443)(m" 1 ) 


Figure 2. Histogram of the total nonwater absorption at 443 nm 
Otw(443), m -1 ) for SBI spr, SBI sum, MR (left) and CASES, NA- 
BOS, MALINA (right). 


ues were calculated using linear regression analysis. To eval- 
uate the performance of our CDOM algorithm, root mean 
square error (RMSE), mean normalized bias (MNB), and ab- 
solute percent difference (APD) were used: 


RMSE = 


[1°Sio^cdom ( 443) log 10 ^Q OM (443)] 

n = 1 


N 


( 1 ) 


2.2 Statistical analyses 

Bio-optical quantities are often log-normally distributed in 
natural environments (Campbell, 1995). The normality of 
distribution for log-transformed a tw ( 443) and <2cdom(443) 
values was examined for Arctic waters using a Kolmogorov- 
Smirnov test. If the normality of distribution was verified 
for a variable, we conducted either a t test (two variables) 
or ANOVA (analysis of variance) to examine a difference in 
mean values for each pair of data. If a variable was not nor- 
mally distributed, we conducted a nonparametric Wilcoxon 
rank-sum test. Geometric mean and geometric SD (standard 
deviation) were thus obtained for these variables in this study. 
Otherwise, arithmetic mean and SD were used, r 2 and p val- 


MNB = x V ^ QglQ<2 CD0M(443) log 10 ^cD QM (443) 

N hi lo glO< D W 443 ) 


£00 ^ A |lOglQfl™pQ M ( 443 ) l°glO^CDOM ( 443 ) | 

N hi lo glO«CDW 443 ) 


3 Results and discussion 
3.1 Absorption budget 

In view of ocean color applications for the Arctic Ocean, 
the total nonwater absorption properties were examined first. 
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Figure 3. Histogram of CDOM absorption at 443 nm Ocdom(443), 
m -1 ) for SBI spr, SBI sum, MR (left) and CASES, NABOS, MA- 
LINA (right). 


We found that atw(443) in river-influenced coastal waters 
(i.e., CASES, NABOS, and MALINA) tended to be higher 
than in oceanic waters (i.e., SBI spr, SBI sum, and MR) 
(Fig. 2). This result is consistent with the fact that river dis- 
charge largely influences absorption properties for the south- 
ern Beaufort and Siberian seas, but not for the Chukchi Sea 
and the western part of the Beaufort Sea (Matsuoka et al., 
2011, 2012; Amon et al., 2012). Over the Siberian shelves 
(or EAO), <2 tw (443) was significantly higher ( p < 0.001) than 
in the southern Beaufort Sea (or WAO) (Fig. 2; Table 1), 
which was also indirectly shown by Aas et al. (2002) using 
spectral diffuse attenuation coefficients of downwelling irra- 
diance (K^k), m -1 ) (the total nonwater absorption was not 
directly measured in their work). 

The higher <2cdom(443) of the EAO was partly responsible 
for the higher a tw ( 443) compared to that of the WAO (rela- 
tive contribution to a tw (443) of 85 zb 7 % and 67 ± 19 %, re- 
spectively; Fig. 3). This is also consistent with the fact that 
the load of DOC and CDOM in the Siberian rivers (e.g., 
Lena, Yenisei, and Ob) is significantly higher than in the 
Mackenzie River (Stedmon et al., 2011). CDOM absorp- 
tion is strongly correlated with DOC concentrations in river- 
influenced coastal waters of the Arctic Ocean (Matsuoka et 



Figure 4. The contribution of absorption coefficients of phyto- 
plankton (a^(443), m -1 ), NAP (a NAP (443), m” 1 ), and CDOM 
("CDOM (443), m -1 ) absorption coefficients to the total nonwater 
absorption at 443 nm. See Fig. 1 for symbols. 


acD0M(443)(m *) 


ao(443){m ^ 


a NA p(443)(m 0 


al., 2012; Walker et al., 2013). The higher CDOM absorp- 
tion in the EAO likely reflects higher DOC concentrations 
in high-latitude watersheds of Siberia compared to those in 
North America (Raymond et al., 2007). This topic is further 
discussed in Sect. 3.4. 

The strong contribution of CDOM absorption is consis- 
tent with the fact that the Arctic Ocean receives the largest 
amount of freshwater relative to its volume (1 1 % of global 
freshwater input while its volume is only 1 % of the global 
ocean, Siklomanov, 1993). While this high contribution of 
CDOM absorption is well acknowledged especially for the 
WAO in both coastal and offshore waters (Belanger et al., 
2006, 2013; Matsuoka et al., 2007, 2009, 2011), it is evi- 
denced as well here for the EAO (Fig. 4), confirming the 
characteristic for the whole Arctic region (68 ±19%, N = 
603; Table 1). Consequently, waters include a high propor- 
tion of CDOM for both oceanic and coastal waters of the Arc- 
tic Ocean, significantly influencing several estimates such as 
K&(k) and PP as derived from ocean color data (e.g., Arrigo 
and van Dijken, 2011; Reynolds et al., 2014; Arrigo et al., 
2014). 

3.2 Parameterization 

While Matsuoka et al. (2011) examined absorption proper- 
ties of the WAO, the majority of the data in that study were 
obtained for oceanic waters. To generalize the absorption 
properties for Arctic waters, large data sets of absorption in 
coastal waters from both the southern Beaufort and Siberian 
seas were added and analyzed in this study. 

To examine differences in bio-optical properties between 
the WAO and EAO, a p ( 443) was plotted as a function of chi a 
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chi a(mg m 3 ) 



Figure 5. (a) Relationship between the total particulate absorption 
coefficients at 443 nm (a p ( 443), m -1 ) and chi a for our Arctic data 
sets. The regression line (d p ( 443) = 0.0548[chl a] 0 658 , N = 657; 
p < 0.0001) and ±95 % confidence intervals are displayed. The re- 
gression line obtained by Bricaud et al. (1998) (grey) is also over- 
laid for comparison, (b) Relationship between <2^(443) and chi a 
(ay (443) = 0.03 15 [chi d] 0J33 ,N = 657 ; p < 0.0001). 


(Fig. 5a). Even when the outliers, observed outside the 95 % 
confidence intervals, were included, a p ( 443) showed signifi- 
cant positive correlation with chi a (; r 2 = 0.62, p < 0.0001; 
N = 657), which was quite similar to the correlation ob- 
served by Bricaud et al. (1998) at lower latitudes. A high 
pigment packaging effect has been acknowledged as a com- 
mon feature for Arctic phytoplankton absorption, resulting 
in lower chi a-specific a<p(k) (a* (A.), m 2 mg chi a~ l ) values 
relative to those at lower latitudes (e.g., Mitchell, 1992; Mat- 
suoka et al., 2011). Because particulate absorption is com- 
posed of phytoplankton and NAP absorption, higher NAP ab- 


Table 2. Average values and standard deviations of «nap( 443) nor- 
malized by SPM («na p(443)*,m 2 g -1 ). For comparison, values ob- 
tained for the COASTLOOC cruises in the literature (Babin et al., 
2003) are also shown. 



Cmise 

«NAP(443)* 

N 

This study 

CASES 

0.0360 ±0.0275 

107 


NABOS 

0.0014 ±0.0019 

15 


MALINA 

0.0602 ±0.0300 

121 


All 

0.0459 ±0.0324 

243 


All but NABOS 

0.0488 ±0.0312 

228 

Babin et al. (2003) 

COASTLOOC 

0.0410 ±0.0230 

328 

Table 3. Correlations between a^ik) and a 0 (443) in this study. 

dfj) (k) at wavelength 412 490 

555 670 

N 

r 2 

0.98 0.99 

0.87 0.95 

993 


sorption values for Arctic waters compared to those at lower 
latitudes (p < 0.05; Table 2) leads to a similar a v ( 443) ver- 
sus chi a relationship. The higher NAP likely originates from 
materials derived from glaciers and delivered by river runoff 
or from coastal erosion by ice (Mitchell and Holm-Hansen, 
1991 ; Hodgkins et al., 2003), by accumulation of organic ma- 
terial in sea ice due to atmospheric deposition of aerosols 
or by ice algae in situ production (Belanger et al., 2013). 
The SPM-specific NAP absorption at 443 nm («nap(443), 
m 2 g -1 ) for the WAO was significantly higher than for the 
EAO (Table 2). This difference likely results from geological 
differences between the two areas (e.g., Amon et al., 2012). 

Figure 5b represents the relationship between rfy(443) and 
chi a. Less variability in this relationship was observed com- 
pared to the a p (443) versus chi a relationship. A strong cor- 
relation between the two variables was found (r 2 = 0.78, 
p < 0.0001; N = 657). As expected, this regression line at 
chi a < 1.0 mg m 3 was significantly below the one obtained 
by Bricaud et al. (1998) at lower latitudes, demonstrating the 
higher pigment packaging effect for Arctic waters. 

It should be stressed that while significantly higher 
^cdom(443) in the EAO was observed compared to in the 
WAO (Figs. 3, 4), the ^(443) versus chi a relationships be- 
tween the EAO (or NABOS data set) and WAO (or the rest 
of the data sets) showed no significant difference or were sta- 
tistically similar (p = 0.35, F test). Despite the geographical 
difference between the WAO and EAO, this result suggests 
that the same relationship can be applied for both. More im- 
portantly, because our semianalytical CDOM absorption al- 
gorithm is based on chi a -specific a^k) estimates (Matsuoka 
et al., 2013), this result highlights that our algorithm can be 
applied not only for the WAO but also for the EAO. 

Matsuoka et al. (2011) and Wang et al. (2005) showed 
that dy (k) in the visible spectral domain is highly correlated 
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Figure 6. Relationship between <2^(443) and a^iX) at 412, 490, 
555, and 670 nm. Coefficients of determination are shown in Table 
3. 


with <2^(443), which can be obtained using ocean color data. 
While a large data set including both oceanic and coastal 
waters was used in this study ( N = 993), a^iX) showed ex- 
cellent correlation with <2^(443) (Fig. 6, Table 3). This kind 
of statistical approach allows obtaining spectral a ^ (A.) values 
using ocean color data, which can be further used for operat- 
ing a spectral primary production model from space. 

Statistics for <2cdom(443) and <2nap(443) and their spectral 
slopes are shown in Table 4. Results are consistent with those 
provided by Matsuoka et al. (2011), although much larger 
data sets obtained from various locations were used in the 
present study ( N = 859, compared to N = 408 in Matsuoka 
et al., 2011), suggesting the ubiquitous features in waters of 
the Arctic Ocean. 

3.3 Evaluation of the CDOM absorption algorithm 

Estimates of DOC concentrations in Arctic river- influenced 
coastal waters using ocean color data are reliable when 
CDOM absorption values are appropriately derived (Mat- 
suoka et al., 2013). We evaluated the performance of a tuned 
semianalytical CDOM absorption algorithm using our Arc- 
tic data sets (Fig. 7); in this study, all data from MALINA, 
ICESCAPE1 and ICESCAPE2 data sets were used for this 
evaluation, whereas Matsuoka et al. (2013) only used data 
from the MALINA data set. It is important to recall that 
this evaluation was performed using independent data sets 
(e.g., MALINA, ICESCAPE1, and ICESCAPE2 data sets) 
that were not used for developing the algorithm (i.e., data 
from SBI spr, SBI sum, and MR cruises were used for 
this development). Results showed that the semianalytical 
CDOM algorithm works reasonably well (Table 5). On av- 



Figure 7. A comparison of estimated CDOM absorption 
0 ^dq M ( 443), m -1 ) and in situ measurements Ocdom(^ 43), 
m -1 ) using data sets (MALINA, ICESCAPE1, and ICESCAPE2) 
that were not used for developing our CDOM algorithm (i.e., data 
from SBI spr, SBI sum, and MR cruises were used for the algo- 
rithm’s development). Phytoplankton absorption spectra for four 
data points above the 50 % error are shown and explained in Fig. A3 
and Appendix A3. 

erage, <2cdom(443) can be derived using our CDOM algo- 
rithm with a mean normalized bias (MNB) of 8.58 % and ab- 
solute percent difference (APD) of 11.72%. Only four data 
points were located above the 50 % error intervals (circles in 
Fig. 7). A possible cause for these values could be the pres- 
ence of mycosporine-like amino acids (MAAs) released by 
algal cells in CDOM (see Appendix A3). Although in situ 
R rs ( X) data in the EAO are currently not available to val- 
idate the algorithm in this area, this result gives confidence 
that a cdom(443) can be derived accurately using ocean color 
data. 

3.4 Ocean color application 

Based on the evaluation of our CDOM algorithm as 
well as bio-optical properties as shown above, we de- 
rived <2cdom(443) from satellite data (^cdom( 443), m -1 ) 
by applying this algorithm to MODIS climatological data 
(Fig. 8). To our knowledge, this is the first synoptic view of 
^cdom(443) for the whole Arctic Ocean using a semianalyt- 
ical approach. Waters with high CDOM absorption are more 
widely distributed in the EAO (e.g., Laptev and Kara seas) 
than in the WAO (e.g., southern Beaufort Sea). For August, 
geometric mean (G. Mean) and geometric standard deviation 
(G. SD) values of ^cdom( 443) i n the southern Beaufort Sea 
(0.059 and 2.363 m -1 , respectively) were similar to our field 
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Table 4. Statistics of «cdom(443) and «nap(443) at 443 nm with their spectral slopes, ScDOM and Snap* respectively. 


Cruise 

^CDOm(443) 
geometric 
mean (SD) 

S CDOM 

arithmetic mean 
±1 SD 

«NAp(443) 
geometric 
mean (SD) 

Snap 

arithmetic mean 
±1 SD 

N 

SBI spr 

0.0534 (2.1231) 

0.0184 ±0.0036 

0.0114(2.4180) 

0.0095 ±0.0022 

80 

SBI sum 

0.0429(1.7890) 

0.0205 ±0.0035 

0.0115(2.7128) 

0.0102 ±0.0015 

111 

MR 

0.0662 (1.6732) 

0.0155 ±0.0031 

0.0062 (2.1135) 

0.0108 ±0.0009 

179 

CASES 

0.0992 (2.6003) 

0.0198 ±0.0027 

0.0187 (4.3476) 

0.0103 ±0.0011 

132 (173) 

NABOS 

0.2463 (1.8366) 

0.0209 ±0.0044 

0.0141 (2.7972) 

0.0165 ±0.0062 

31 (22) 

MALINA 

0.0559 (2.1947) 

0.0185 ±0.0013 

0.0078 (4.8963) 

0.0070 ±0.0014 

326 (436) 

All 

0.0642 (2.2577) 

0.0184 ±0.0032 

0.0095 (3.9657) 

0.0090 ±0.0026 

859 (1001) 



a cooM(443)(m b 
0.3 


Figure 8. Climatology of CDOM absorption at 443 nm from satellite sensors (^qdom^^), m -1 ) from June to September, obtained by 
applying our CDOM absorption algorithm to MODIS monthly-averaged climatology of R r sW data over the 2002-2012 period. Three areas 
were defined for examination of ^CDOM^^) values as follows: southern Beaufort (SB: 125-145° W, 68-72° N), Laptev (LP: 100-145° E, 
69-77° N), and Kara seas (KR: 50-100° E, 64-77° N). Histograms of ^cdom( 443) values within these areas are shown in Fig. 9. 


Table 5. Comparison of «cdom(443) estimates (a^Q M (443), 
m -1 ) with in situ «cdom(443) (a£^Q M (443), m -1 ) using indepen- 
dent data sets that were not used for developing the CDOM algo- 
rithm. Data from MALINA, ICESCAPE1, and ICESCAPE2 were 
used. RMSE, MNB, APD,r 2 , intercept, and slope are provided. 


Data sets 

r 2 

Intercept 

Slope 

RMSE 

MNB 

APD 

N 

This study 

0.87 

-0.022 

0.97 

0.069 

8.58 

11.72 

79 


observations for the MALINA cruise (0.055 and 2.265 m -1 , 
respectively) during the same time period (Matsuoka et al., 
2012). G. Mean and G. SD values for ^cdom(443) 


Laptev (0.138 and 2.833 m -1 , respectively) and Kara seas 
(0.073 and 3.037 m -1 , respectively) are higher than those in 
the southern Beaufort Sea (Fig. 9). The highest ^cdom( 443) 
value observed in the Kara Sea (4.353 m -1 ) is consistent 
with in situ measurements (up to 5m -1 ; Stedmon et al., 
2011; Heim et al., 2013; Orek et al., 2013). In addition to 
our evaluation shown in Fig. 7, this result highlights the reli- 
ability of satellite ocean color data for Arctic waters (see also 
Fig. 5 in Matsuoka et al., 2013). 

To estimate DOC concentrations in river-influenced 
coastal waters using ocean color data, a DOC concentra- 
tion versus CDOM absorption relationship needs to be es- 
tablished. Matsuoka et al. (2012) recently showed a strong 
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Figure 9. Histogram of ^cdom( 443) values in the SB, the LP, and 
the KR seas defined in Fig. 8. 


correlation between the two variables in southern Beaufort 
Seawaters (r 2 = 0.97; p < 0.0001). Similarly, DOC concen- 
trations were well correlated with <2cdom(443) in the EAO 
(Walker et al., 2013). It is important to note that the DOC 
versus <2cdom(443) relationships differ between the WAO 
and EAO (Fig. 10), suggesting that different relationships 
between DOC and <2cdom(443) need to be applied for the 
EAO and WAO. Different DOC versus salinity relationships 
between the WAO and EAO were also reported (Hansell et 
al., 2004). There are likely differences in vegetation between 
watersheds that influence the relative amounts of colored and 
noncolored dissolved organic matter. 

To estimate DOC concentrations in the near surface wa- 
ters of the Arctic Ocean, the regression of DOC versus 
^cdom(443) shown in Fig. 10 was applied to the 4 dom( 443) 
image. Similarly to ^cdom( 443) distribution, waters show- 
ing high DOC concentrations were widely distributed in the 
EAO compared to WAO (Fig. 11). The G. Means of DOC 
concentrations estimated from satellite data (DOC sat , pM) in 
the southern Beaufort Sea (82 pM) were of the same order of 
magnitude as in situ measurements for the MALINA cruise 
(79 pM: Matsuoka et al., 2012). The values in the Laptev 
(279 pM) and Kara seas (267 pM) were higher than in the 
southern Beaufort Sea (Fig. 12), which is consistent with in 
situ measurements (Amon et al., 2012; Walker et al., 2013). 
The high DOC concentrations in the EAO are likely due to 
higher DOC concentrations in the high-latitude watersheds 
of Siberia compared to those in North America (Raymond et 
al., 2007; Stedmon et al., 2011). 

As an example of application of our DOC algorithm, spa- 
tial variability in DOC estimates from 2009 to 20 1 1 was ex- 
amined (Fig. 13). While the Mackenzie River plume stayed 
near coast on 4 August 2010, the plume extended far offshore 
toward north approximately two weeks later. This plume di- 
minished on 7 September 2010. In 201 1, river plume tended 



Figure 10. In situ relationships between DOC and «cdom(443) 
for the WAO (M12: Matsuoka et al., 2012; r 2 = 0.97, p < 
0.0001; DOC = 55 + 357 • a C DOM(443)) and EAO (W13: 
Walker et al., 2013; r 2 = 0.90, p < 0.0001; DOC = 245 + 
171 • «cdom(443)). Data points and 95% confidence intervals 
for M12 are shown as black circles and red dotted lines, re- 
spectively. The ranges of <2 cdom( 443) and DOC concentrations 
for the regressions provided by Ml 2 and W13 and applied in 
this study were as follows: 0.018 m -1 <«cdom(443)< 1.08 m -1 ; 
55 pM<DOC<500 pM; 0.39 m -1 <«cdom(443)< 8.4 m -1 ; and 
166 pM < DOC < 1660 pM. 


to distribute North West direction from 29 July to 2 August 
2011. The plume then diminished by 12 September 2011. 
Note that patterns of the Mackenzie River plume distribu- 
tion are mainly controlled by wind direction (e.g., Carmack 
et al., 1989). Thus, spatial and temporal variability in DOC 
estimates were clearly observed using our algorithm. 

Ongoing global warming will likely increase the release 
of presently sequestered carbon that originates, in part, from 
thawing of the permafrost (e.g., Frey and Smith, 2005; Sted- 
mon et al., 2011). Continuous monitoring and quantification 
of DOC concentrations and the budget for the whole Arctic 
Ocean are therefore urgently required to better understand 
modifications in carbon cycling as a result of global warm- 
ing. 


4 Conclusions 

This study examined light absorption properties of the Arc- 
tic Ocean using a large data set including both oceanic and 
coastal waters. Our results clearly showed the specificity of 
those properties: high proportion of CDOM absorption, and 
high pigment packaging effect. Despite the high proportion 
of CDOM absorption to the total nonwater absorption, the re- 
lationship between phytoplankton absorption at 443 nm and 
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Figure 11. Climatology of DOC concentrations in Arctic river- influenced coastal waters (DOC sat ) from June to September, obtained by 
applying the DOC versus «cdom(443) regressions to the ^cdom( 443) images shown in Fig. 8. For the WAO, the regression obtained by 
Matsuoka et al. (2012) was applied within the measured range of «cdom(443) (i.e., 0.018<acdom(443)< 1.08 m -1 ). For the EAO, the 
regression obtained by Walker et al. (2013) was applied within the measured range of «cdom(443) (i.e., «cdom(443) < 8.4 m -1 ). To obtain 
those values, «cdom(443) was estimated from «cdom(350) by assuming ScDOM — 0.0175 nm _1 for the EAO (Aas et al., 2002; Stedmon 
et al., 2011). Histograms of DOC sat values within the white boxes defined in Fig. 8 are shown in Fig. 12. 
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Figure 12. Histogram of DOC sat values in the SB, LP, and KR seas 
defined in Fig. 8. 


chi a concentration was not significantly different between 
the WAO and EAO. Our semianalytical CDOM algorithm 
depends on chi ^-specific a^iX) estimates, highlighting the 
applicability for the whole Arctic Ocean. Using DOC ver- 
sus < 2 cdom( 443 ) relationships, DOC concentrations in river- 
influenced coastal waters can now be estimated semianalyt- 


ically using satellite ocean color data. Our statistical results 
demonstrated that both <2cdom(443) values and DOC con- 
centrations obtained using satellite ocean color data are rea- 
sonable compared to in situ measurements, suggesting its po- 
tential utility for obtaining a quantitative estimate of the car- 
bon budget of the Arctic Ocean. 

The DOC budget of the Arctic Ocean is expected to 
change rapidly due to sea ice reduction, permafrost thaw- 
ing, and river discharge increases as consequences of global 
warming. To examine this issue temporally and geographi- 
cally using satellite ocean color data, we developed a semi- 
analytical algorithm for estimating DOC concentrations in 
Arctic river-influenced coastal waters. Compared to empiri- 
cal algorithms that have been proposed (e.g., Mannino et al., 
2008; Fichot and Benner, 2011), because CDOM absorption 
is semianalytically derived from ocean color data, this algo- 
rithm depends less on empirical relationships established for 
particular time periods and areas. Another semianalytical al- 
gorithm using difference equations for inherent optical prop- 
erties (IOPs) was also developed in the Kara Sea (Korosov et 
al., 2012). Our pan- Arctic data set and robust semianalytical 
methods should improve satellite estimates of DOC concen- 
trations in surface waters of the Arctic Ocean and help to 
enhance our understanding of modifications of DOC cycles 
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Figure 13. Satellite-derived DOC concentrations in the surface layer for selected MODIS Aqua satellite data recorded in 2009 (top), 2010 
(middle), and 2011 (bottom). 


in the Arctic Ocean. Further work regarding interannual vari- 
ability in DOC budget is required. 

We acknowledge that an issue regarding the seasonal vari- 
ability in the relationship between DOC concentrations and 
CDOM absorption needs to be considered. While no study 
examining this issue has yet been conducted so far for Arc- 
tic waters, Mannino et al. (2008) demonstrated in the U.S. 
middle- Atlantic Bight that the intercept of the DOC versus 
^cdom(443) relationship changed seasonally in their study 
region, whereas the slope remained unchanged. This result 
suggests that seasonal relationships should be established to 
improve the estimates of DOC concentrations from space. 
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Appendix A: 

Al Comparison of chi a determination by the HPLC 
and fluorometric methods 

Chi a concentrations determined fluorometrically (chi a fluo ) 
were compared to those using the HPLC (chi a R?LC ) tech- 
nique (Fig. Al). Correlation between the two quantities was 
high (r 2 = 0.89, N = 177), and slope was not significantly 
different from the 1:1 line ( p < 0.0001). This result sug- 
gests that in this study, chi < 2 fluo can appropriately replace 
chi a HPLC when chi a HPLC is not available. 



Figure Al. Comparison of chi a determined fluorometrically (chi 
cfl uo ) and chi a using the HPLC method (chi a HPLC ). 


A2 Evaluation of estimated phytoplankton absorption 
using the Bricaud and Stramski (1990) method 

We evaluated the performance of the Bricaud and Stramski 
(1990) method to estimate phytoplankton absorption using 
total particulate absorption data from the ICESCAPE1 cruise 
(Fig. A2). Results showed that a<p{X) values can be derived 
with high accuracy using this method at wavelengths in the 
visible spectral domain. This method was used for the NA- 
BOS data set. 



Figure A2. Comparison of a^iX) estimated using the Bricaud and 
Stramski (1990) method with in situ measurements during the 
ICESCAPE1 cruise. 


A3 High UV absorption for some phytoplankton 
communities 

High values in phytoplankton absorption spectra around 
330 nm were observed at four stations during ICESCAPE 
cruises. The presence of mycosporine-like amino acids 
(MAAs) in natural waters has been reported at all latitudes 
(e.g., Karsten et al., 1991, 1998; Whitehead and Vemet, 
2000; Tilstone et al., 2010). Most of the MAAs’ absorp- 
tion peaks are observed at 330-334 nm due to several types 
of MAAs (e.g., shinorine, porphyra-334, asterina; Carreto 
et al., 2005; Llewellyn and Airs, 2010). A clear bump was 
also observed at 360 nm at station 38 during ICESCAPE 1, 
which corresponded to the peak of palythene (Llewellyn and 
Airs, 2010). Although actual concentrations of MAAs dur- 
ing our observations are not available, the strong UV (ultra- 
violet) absorption observed for some samples in this study 
suggest MAAs were present. Since MAAs are water sol- 
uble, they may be released from the cells upon freezing- 
melting processes of sea ice, senescence, or grazing enter- 
ing the dissolved organic pool and hence contribute to our 
<2cdom(443) estimates (circles in Fig. 7). It is interesting to 
note that a^iX) at 330 and 360 nm can influence R rs ( X) in 
the visible spectral domain. A similar result was also ob- 
served using an empirical relationship between <2cdom(443) 
and /Ui(380)//fd(780) (Hooker, unpublished). The strong im- 
pact of MAAs’ absorption and spectral reflectance in the 
UV spectral domain during a massive red tide was reported 
by Kahru and Mitchell (1998) and is the basis for the red 
tide index published for the global imager (GLI) algorithm 
(Mitchell and Kahru, 2009). This influence of MAAs on 
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UV-visible absorption may allow estimates of MAAs using 
ocean color data in the future if the challenge of atmospheric 
correction in the UV spectral domain can be solved. 



Figure A3. Spectra of phytoplankton absorption coefficients 
(< d(p(X . ), m -1 ) at four stations during ICESCAPE1 and 2 cruises. 
CDOM absorption estimates at these points were outside of the 
50 % errors as shown in Fig. 7. 
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